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ABSTRACT 

We investigate the formation and evolution of circumstellar disks in turbulent cloud 
cores until several 10^ years after protostar formation using smoothed particle hydro- 
dynamics (SPH) calculations. The formation and evolution process of circumstellar 
disk in turbulent cloud cores differs substantially from that in rigidly rotating cloud 
cores. In turbulent cloud cores, a filamentary structure appears before the protostar 
formation and the protostar forms in the filament. If the turbulence is initially suffi- 
ciently strong, the remaining filament twists around the protostar and directly becomes 
a rotation-supported disk. Upon formation, the disk orientation is generally misaligned 
with the angular momentum of its host cloud core and it dynamically varies during 
the main accretion phase, even though the turbulence is weak. This is because the 
angular momentum of the entire cloud core is mainly determined by the large scale 
velocity field whose wavelength is comparable to the cloud scale, whereas the angular 
momentum of the disk is determined by the local velocity field where the protostar 
forms and these two velocity fields do not correlate with each other. In the case of 
disk evolution in a binary or multiple stars, the disks are misaligned with each other 
at least during the main accretion phase, because there is no correlation between the 
velocity fields around the position where each protostar forms. In addition, each disk 
is also misaligned with the binary orbital plane. Such misalignment can explain the 
recent observations of misaligned disks and misaligned protostellar outflows. 

Key words: star formation - circumstellar disk - protoplanetary disk - planet forma- 
tion — binary system - methods: hydrodynamics - smoothed particle hydrodynamics 



1 INTRODUCTION 

Many circumstellar disks are observed in star-forming re- 
gions. They are by-products of star formation and are di- 
rectly connected to planet formation. Thus, it is important 
for understanding star and planet formation to clarify the 
formation and evolution of circumstellar disks. Recent stud- 
ies on protostar and disk formation in a collapsing molecu- 
lar cloud core suggest that the circumstellar disk forms at a 
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ver y early phase of star formation (Bate 1998; Walch et al. 
|m)9; Machida et al. 2010; Inutsuka et al. 2010; Bate 2011). 
and is more massive than the protostar for at lea st 10^ 
years after protostar formation (Inutsuka et al.ir201Ql V Dur- 
ing this phase, the circumstellar disk is gravitationally un- 
stable and develops non-axisymmetric spiral arms. However, 
these studies have examined the formation and evolution of 
the circumstellar disk in the rigidly rotating cloud core. 

Molecular clouds, which typically have a scale > 0.1 pc, 
usually display complex internal motions that are observed 
as a broad emission-line profile. The internal motion in a 



molecular cloud is believed to be caused by turbulence. Al- 
though a detailed generation process of the internal motions 
is still unclear, observed line profiles in molecular clouds are 
consistent with Gaussian veloc i ty field s with a Kolmogorov 
spectrum dPubinski et al ] |l995l : lKlessenll20 00). Observations 
also indicate the existence of turbulence even for molecu- 
lar cloud cores, which typically have a scale of 0.1-0.01 pc. 
iGoodman et al.l (|l99 3') showed that the specific angular mo- 
mentum of a molecular cloud core is roughly proportional 
to 2 i?^"^, which is derived from observations of molec- 
ular cloud cores with a size of 0.06 — 0.6 pc. This scaling 
law is in good agreement with the scaling, j oc i?^"^ which 
is produced by a turbulence field with velocity power spec- 
trum, P(k) oc . Note that when the molecular cloud core 
rotates rigidly, the specific angular momentum should obey 
the scaling law of j oc . iBurkert h Bodenheimerl (|2000l ) 
pointed out that random Gaussian velocity fields with power 
spectra of P(A:) oc can reproduce the observed projected 
rotational properties of molecular cloud cores. 
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The evolution of a circumstellar disk during Class 0- 
I phases may be strongly affected by the velocity field of 
the molecular cloud core because it is mainly determined 
by mass accretion from the (turbulent) envelope. Thus, it 
is important to investigate the disk evolution in turbulent 
molecular cloud cores for a comprehensive understanding of 
the evolution of circumstellar disks. 

In addition to formation and evolution of the circum- 
stellar disk around a single protostar, cloud turbulence is 
also expected to affect the formation and evolution of cir- 
cumstellar disks in binary systems. In a collapsing cloud 
core, fragmentation may occur, causing a binary system to 
appear. In a rigidly rotating cloud core, disks in the binary 
system are aligned with each other by definition and also 
with the binary orbital plane (see, iTsukamoto Machidal 
l201lh . On the other hand, observational evidenc e indi- 
cates m isalignment of binary disks. For example, Dav is et al.l 
(|l994l l observed pr otostellar jets with d i fferen t orientations 
in a binary system. iRoccatagliata et al ] (|201lh recently ob- 
served two disks in a proto-binary system. In their obser- 
vation, the primary protostar has an almost face-on disk, 
whereas the secondary disk is edge-on toward the observer. 
Thus, the disk orientation around the primary protostar is 
almost perpendicular to the orientation of the disk around 
the secondary protostar. With observations of many wide 
binary systems, ' Hald ([l994) showed that the stellar rota- 
tional equatorial plane is often misaligned with the binary 
orbital plane. These observations indicate that in wide bi- 
nary systems, the disk plane is frequently misaligned with 
the binary orbital plane. However, simulations starting from 
a molecular cloud core with simple systematic rotation can- 
not explain these observations. 

The gravitational collapse of a turbulent molecu- 
lar cloud c ore has been studied by several groups. 
iMatsumoto Hanaw3 (|201ll ) calculated the gravitational 
contraction of turbulent cloud cores and showed the mor- 
phological evolution of the collapsing cloud core before 
protostar formation. With long term cal culations (until 
~ 10^ years after the protostar formation), I Goodwin et al.l 
(|20Q4 b.a) showed that fragmentation frequentl y occurs even 
i n a w eakly turbulent environment. Recently, IWalch et al.l 
(|2Q1QI 1 calculated the evolution of turbulent cloud cores and 
showed that protostar and circumstellar disk formation in 
turbulent cloud cores is different from that in rigidly rotat- 
ing cloud cores. They showed that protostars form from a 
filamentary structure and pointed out that the formation 
condition for a binary or multiple stellar system rarely de- 
pend s on the tot al angular momentum of the host cloud 
core. IWalch et al.l (|2012i ) investigated the evolution of low- 
mass cold cores. They focused on the relationship between 
the evolution of the core and the maximum wavelength (or 
minimum wavenumber) of the turbulence and found that 
the maximum wavelength of the turbulence is a critical pa- 
rameter for the evolution of cloud cores. In their simulation, 
the turbulent energy is fixed and the dependence of the disk 
evolution on the turbulent energy is still unclear. 

In this study, we investigate circumstellar disk for- 
mation in a turbulent cloud core over > 10^ years. 
We have already reported the disk formation in cloud 
cores with systematic rotation (or rigid rotation) in 
ITsukamoto Machidal (|201ll ) in which we calculated the 
cloud evolution with different classical cloud parameters 



(e.g., iMivama et al.l [l984l ) representing the cloud thermal 
(a) and rotational (/3) energ ies. In the present study , whic h 
is complementary to that of 'Tsukamoto & Machidaj (|201lh . 
we investigate the effects of turbulence strength on the evo- 
lution of the protostar and circumstellar disk with one cloud 
parameter representing thermal (a) energy and one repre- 
senting turbulent (7turb) energy. We also investigate disk 
evolution and its orientation in a binary system. This pa- 
per is organized as follows. In §2, we describe the numerical 
method and initial conditions. In §3, we present the numer- 
ical results. Finally, in §4, we discuss our results. 



2 NUMERICAL METHOD AND INITIAL 
CONDITIONS 

2.1 Numerical Method 

Our simulations were conducted using the smoothed particle 
hydrodynamics (SPH) code, which we used in our previous 
study (|Tsukamoto Machidal 1201 li V The code includes an 
individual time-step technique and uses the Barnes-Hut tree 
algorithm to calculate the self-gravity with an opening an- 
gle = 0.5. We include an artificial viscosity according as 
prescribed by Monaghan (1997) with = 1 and also use 
the Balsara switch (|Balsarall 19951 ) . Our code was parallelized 
with MPI and was verified using several standard test prob- 
lems. 

To mimic the t hermal ev olution of the cloud core calcu- 
lated by Masunag a Inutsu ka (2000), we adopted the fol- 
lowing barotropic equation of state. 
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where Cs,o — 190 m s~^ and pc 

In addition, to calculate the evolution of the circumstel- 
lar disk for ^10^ ^years, we adopted the sink particle tech- 
nique described by Bate et al. (1995"). We assume that a pro- 
tostar forms when the particle density exceeds the threshold 
density, psink = 4 x 10~^g cm~^. Next, a sink particle with 
an accretion radius of 1 AU is dynamically introduced. 



2.2 Initial Settings 

As the initial state, we adopt a spherically-symmetric cloud 
core with an isothermal temperature of T = 10 K. Each 
cloud core has a uniform density within the range pinit — 
6.9 X 10~^^ to 1.9 X 10~^^cm~^ and a size ranging from 
R = 1967 to 5900 AU (see. Table 1). Ah models have the 
same cloud mass of IMq. The initial cloud cores are modeled 
with about 520000 SPH particles. The mass resolution of all 
calculations is 1.9x 10~^ Mq and our calculations fulfill the 
resolution requirement suggested by Bate & Burkert (19971). 
The smoothing length of i th particle, hi is given as 



hi 1.2 



rrii 



1/3 



(2) 



where m^, rhoi are the mass and the density of i th parti- 
cle. Only a turbulent velocity field without the systematic 
rotation velocity is imposed on the initial cloud core. 
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Table 1. Model parameters 



Model 


OL 


7turb 


R (AU) 


/^eff 


Pinit (g cm 3) 


Mach number 1 


Fragmentation 


1 


0.6 


0.1 


5900 


1.3 X 10-2 


6.9 X 10-19 


0.66 


N 


2 


0.6 


0.06 


5900 


7.9 X 10-3 


6.9 X 10-19 


0.50 


N 


3 


0.6 


0.03 


5900 


4.1 X 10-3 


6.9 X 10-19 


0.36 


N 


4 


0.4 


0.3 


3933 


4.1 X 10-2 


2.3 X 10-18 


1.4 


N 


5 


0.4 


0.1 


3933 


1.4 X 10-2 


2.3 X 10-is 


0.80 


N 


6 


0.4 


0.06 


3933 


8.5 X 10-3 


2.3 X 10-is 


0.63 


N 


7 


0.2 


0.3 


1967 


4.0 X 10-2 


1.9 X 10-1'^ 


2.0 


Y 


8 


0.2 


0.1 


1967 


1.3 X 10-2 


1.9 X 10-1'^ 


1.1 


Y 


9 


0.2 


0.06 


1967 


8.2 X 10-3 


1.9 X 10-1'^ 


0.87 


Y 



Mach number is the root mean square of the turbulent velocity in the initial cloud core. 



2.2.1 Model Parameters 

Gas accretion from the infalling envelope onto the circum- 
stellar disk controls the evolution of the circumstellar disk. 
The gas accretion rate is related to the thermal and kinetic 
(or turbulent) energies of the host cloud core. In this study, 
to investigate disk evolution in cloud cores with different 
thermal and turbulent energies, we use two parameters, a 
and 7turb- The parameter a is the ratio of thermal energy 
(^thermal) to gravitational energy (^grav) in the initial cloud 
core and given as 



-^'thermal S-RoCg 



l^gravl 2GM 



(3) 



where Rq and M are the initial radius and mass of the cloud 
core, respectively. The parameter 7turb is the ratio of turbu- 
lent energy (^turb) to gravitational energy (^grav): 

-E-turb 



7turb ■ 



(4) 



In equation Q, ^turb and ^grav are numerically estimated 
in the initial cloud core. The strength of the turbulent ve- 
locity field in the initial cloud core can be described by the 
parameter 7turb. 



2.2.2 Turbulence Realization 



iBurkert k, Bodenheimerl (|2000l ) showed that random Gaus- 
sian velocity fields with power spectra of P(k) oc to 
can reproduce the observed projected rotational properties 
of molecular cloud cores. Thus, a molecular cloud core ap- 
parently has a systematic rotation even when it has only a 
turbulent velocity field. We adopt a turbulent velocity field 
with power spectra of P{k) oc A:~^, in which we assume 
a divergence- free velocity field. The prescription for realiz- 
ing turbulence in our initial settings is as follows: First, we 
generate a random Gaussian field with the power spectrum 
Pa (A:) =< \Ak\^ >oc . Next, we compute the Fourier 
transform of the velocity field Vh as 



Vh ^ ik X Ak. 



(5) 



Finally, the velocity field is generated on a 128'^ uniform 
grid, and the particle velocities are interpolated from the 
grid. We adopted the minimum wave number to be /cmin = 1 
and the maximum wave number to be kmax — 128. Unlike 
a rigidly rotating cloud core, the turbulent velocity field is 



not uniquely identified by the parameter set, (a, 7turb) be- 
cause of its stochastic nature. Thus, it is difficult to sys- 
tematically investigate the evolutio n of the turbule n t cloud 
core with parani eters 7turb and a. iGoodwin et alj (|2004al l 
and lWalch eraD (|2012l ) resolved this difficulty by consider- 
ing the ensemble of simulations. Although their approach 
may be advantageous for investigating the turbulent cloud 
core, it requires too many computational resources. There- 
fore, we adopted a different approach for this study in which 
we simulate the evolution of the turbulent cloud core having 
a "typical" angular momentum for each parameter set. To 
construct a velocity field with typical angular momentum 
for a given parameter set (a, 7turb), we generated 1000 dif- 
ferent velocity fields for each parameter set with different 
random number seeds. Next, from the 1000 velocity fields, 
we selected the velocity field whose angular momentum was 
the closest to the mean value of the angular momentum 
and used it to calculate the cloud evolution. This procedure 
makes it possible to simulate a turbulent cloud core that 
has a plausible angular momentum for a given parameter 
(tt, 7turb)- To relate the turbulent energy to the rotational 
energy of the initial cloud core, we introduce the parameter 
/3eff as 



25 J^loud 

12 GM^R' 



(6) 



where Jcioud is the total angular momentum of the cloud 
core. This parameter can be regarded as the effective ra- 
tio of rotational energy to gravitational energy of the cloud 
core. The values of /3eff for each model are listed in Table 1. 
According to this definition, the effective rotational energy 
is roughly related to the turbulent energy as /3eff ^ 0.l7turb- 
Using the parameters a and 7turb, we constructed 9 
models and calculated the evolution of the circumstellar disk 
for each model. The model names and parameters are listed 
in Table 1. 

Note that we adopted the turbulent velocity field but 
the uniform density field as the initial state. Thus, the ve- 
locity field and density field of our initial conditions are not 
self-consistent. We will discuss how this inconsistency would 
affect our results in ^4.5.11 

Note also that some models have somewhat unrealistic 
parameter. Model 7 have large initial Mach number that is 
much larger than typically observed and model 8 and 9 is 
highly gravitationally unstable (a + 7turb = 0.26,0.3). We 
will discuss the validity of these initial conditions in ^4.5.21 
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Figure 1. Time sequence of the logarithm of the face-on surface density before and after protostar formation for model 4 (a = 0.4 and 
7turb = 0.3). The 2;-axis (i.e., the line-of-sight direction) is set parallel to the angular momentum of the entire initial cloud core. The top- 
left and -middle panels show snapshots approximately 2.6 x 10^ and 1.8 x 10^ years before protostar formation, respectively. The top-right 
panel shows a snapshot just when the protostar forms. The bottom left, -middle, and -right panels show snapshots 1.1 x lO'^, 2.1 x 10^, 
and 4.1 X 10^ years after protostar formation, respectively. The elapsed time in the calculation is shown in each panel. 



3 RESULTS 

In this section, we first give an overview of evolution process 
of circumstellar disk in turbulent core and how it changes 
in response to the strength of turbulence. We adopt model 

4 (a = 0.4 and 7turb = 0.3) as a typical model with strong 
turbulence and model 6 (a = 0.4 and 7turb = 0.06) as a 
typical model with weak turbulence, and investigate them 
both in detail. It is observed that whether the initial velocity 
field is supersonic (M > 1) or subsonic (M < 1) is not im- 
portant for the evolution processes. Next, we investigate the 
evolution of the disk orientation for all models and, finally, 
we investigate the formation process of binary or multiple 
systems in a turbulent cloud core and the evolution process 
of circumstellar disks in a binary system. 

3.1 Overview 

3.1.1 Disk evolution in strong turbulent core 

In this subsection, we show how the circumstellar disk 
evolves in a turbulent cloud core that initially has relatively 
strong turbulence. As a typical case, we chose model 4 with 
parameters a = 0.4 and 7turb = 0.3 and an initial root mean 
square (rms) Mach number of 1.4. Figure [T] shows the time 
evolution of the surface density around the center of the 
cloud core. In the figure, the 2;-axis is chosen to be parallel 
to the angular momentum of the entire initial cloud core. 

Before pr o tostar formation, the first ( adiabatic) core 
([Larsonl Il969l : iMasunaga Inutsukal |200Q| ) appears. The 
central region of high-density gas in the top-left and -middle 



panels of Figure [T] corresponds to the first core. These pan- 
els show that filamentary structure accompanies the first 
core, which then, coils around the first core and changes di- 
rectly to a circumstellar disk. The second collapse occurs 
at t = 75500 years (top-right panel). By this epoch, the 
two spiral arms already develop. Although we cannot rec- 
ognize the filamentary structure in the top-right panel, a 
weak large-scale filament remains around the disk at this 
epoch. In addition, the gas accretion from the filamentary 
structure onto the disk continues. After protostar forma- 
tion, this large-scale filament (the structure which resides at 
^ 200 AU) coils around the disk and directly connects to the 
disk's spiral arm (bottom- left panel). Finally, as seen in the 
bottom-right panel of Figure [U the filamentary structure 
accretes onto the disk and disappears. 

Previous studies, which assumed the rigidly rotating 
cloud core with large rotational energy (/3 of the rigidly ro- 
tating core is comparable to /3eff of the turbulent core) , have 
shown that the first core has enough angular momentum to 
develop a bar-mode instability and trailing spiral arms ap- 
pears. The spiral arms remove the angular momentum from 
the central region. With the angular momentum transfer, 
the central density becomes high and the second collapse 
occurs. Then, the remnant of the first core which increases 
in siz e with the spiral arms becomes a circums tellar disk. 
{ Batei Il998l . l201ll : iMachida Matsumotol l201lh . However, 
the disk formation mechanism seen in Figure [1] differs from 
that found in the previous studies. 

The first core formed in the turbulent cloud core has 
less angular momentum than that formed in the rigidly ro- 
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Figure 2. Radial distribution of surface density (top) and radial velocity (bottom) for model 4. Lines corresponds to 1.8 x 10"^ years 
(solid line) before protostar formation, years (dashed line), 1.1 x 10^ years (dotted line), and 4.1 x IC^ years (dashed-dotted line) after 
protostar formation. 



tating cloud core because the gravitational collapse tends 
to occur at a stagnation point of the velocity field around 
which the rotational energy of the gas is relatively low. As 
a result, the first core cannot develop the bar- mode insta- 
bility even with the strong turbulence. Instead, the filamen- 
tary structure, which forms before the first core formation, 
twists around the first core or protostar and, then, becomes 
directly a rotationally supported disk. 

Figure [2] shows the radial profiles of the surface density 
and the radial velocity, for which each value is azimuthally 
averaged. The origin is set at the position of the protostar 
(or the position of the gas particle that has maximum den- 
sity before the second collapse). The velocity is measured 
in the standard of rest of the protostar. The solid line in 
the top panel shows the surface density profile ^ 1.8 x 10"^ 



years before the second collapse, which corresponds to the 
top-middle panel of Figure [H At this epoch, the disk sur- 
rounding the first core ranges from 10 to 50 AU. As described 
above, the filamentary structure, which formed before pro- 
tostar formation, changes directly into the disk. The dashed 
line in the top panel of Figure [2] is the surface density pro- 
file at which the second collapse just begins. The disk radius 
reaches ^ 100 AU at this epoch. The dotted line shows the 
surface density profile at ^ 1.1 x 10^ years after the second 
collapse. At this epoch, as seen in the bottom panel of Fig- 
ure [2 a positive radial velocity arises at approximately 200 
to 300 AU. This flow is caused by the strong m — 1 mode 
of the spiral arm which is generated approximately 100 AU. 
This type of outward flow appears repeatedly and readjust 
the surface density toward the stable configuration. 
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Figure 3. Evolution of Fourier amplitudes of surface density at 100 AU against elapsed time after protostar formation for model 4. 



As shown in Figured] the structure of the disk is highly 
non-axisymmetric, which is caused by non-axisymmetric ac- 
cretion from large-scale filaments. To clarify this, we plot 
in Figure [3] the Fourier amplitudes of the surface density at 
100 AU for 2.0 x 10^ years after protostar formation. The 
Fourier amplitudes are calculated from a Fourier series, 



- Am exp[-im(/)], (7) 



where m is the azimuthal wavenumber and Am is the Fourier 
amplitude. Figure [3] shows that most of the time, the m = 1 
mode dominates the other modes, although the m — 2 mode 
occasionally dominates the m = 1 mode. In addition, the 
amplitude of other higher modes (m > 3) is smaller than 
that of modes m — 1 and 2. The prominent m — 1 mode is 
a significant feature of disks with non-axisymmetric envelope 
accretion. 

Figure |4] shows the time evolution of an edge-on view of 
the center of the cloud core for model 4. The figure indicates 
that the disk orientation varies during main accretion phase. 
We discuss the evolution of the disk orientation further in 

E21 

As we showed above, the evolution process of a circum- 
stellar disk in a strong turbulent core is different from that 
in a rapidly rotating core in many respects. As we will see 
below, the evolution process becomes similar to that of a 
rigidly rotating cloud core as the turbulent energy decreases. 
However, note that the dynamical change of the disk orien- 
tation also occurs even in the weak turbulent core. 



3.1.2 Disk evolution in weak turbulent core 

As a typical case for the formation and evolution of a cir- 
cumstellar disk in a weakly turbulent cloud core, we chose 
model 6 that has parameters oi a — 0.4 and 7turb = 0.06. 
The initial rms Mach number is 0.63. Compared with model 
4 (strong turbulent case), model 6 has the same thermal 
energy but a smaller turbulent energy. 

Figure [5] shows a face-on view of the time evolution of 
the center of the cloud core. Before protostar formation (top- 
left panel), some filaments appear and form a complicated 
structure. However, the density amplitude of filaments ap- 
pearing in this model is much weaker than what appears in 
the models which initially have strong turbulence. As the 
mass accretion onto disk proceeds, the circumstellar disk 
becomes massive and develops the spiral arms (bottom pan- 
els). Finally, the radius of the circumstellar disk reaches 
> 100 AU. The filaments that appear in the early evolution 
phase (top panels) are attributed to anisotropic gas accre- 
tion in a weakly turbulent environment, whereas the spirals 
that appear in the later evolution phase (bottom panels) 
develop because of gravitational instability of the disk itself. 

Figure [6] shows the disk evolution when viewed edge- 
on. The results indicate that the disk orientation varies with 
time even with the weak turbulence. As we will see in ^3.21 
the dynamical change of the disk orientation generally oc- 
curs independent of the strength of the turbulence. The fluc- 
tuation of the disk orientation is a significant feature of disk 
evolution in turbulent cloud cores. 

Figure [71 shows the radial profiles of the surface density 
(top) and radial velocity (bottom) for model 6; the solid, 
dashed, dotted, and dashed-dotted lines correspond to the 
top-middle, bottom- left, bottom-middle, and bottom-right 
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Figure 5. Time sequence of logarithm of face-on surface density before and after protostar formation for model 6 (o; = 0.4 and 7turb = 
0.06). The 2;-axis is parallel to the angular momentum of the entire initial cloud core. Top-left panel shows a snapshot approximately 
4.3 X 10^ years before protostar formation. Top-middle panel shows a snapshot just when the protostar forms. Top-right, bottom-left, 
-middle, and -right show snapshots 8.7 x 10^ years, 3.9 x 10^ years, 1.1 x 10^ years, and 2.6 x 10^ years after protostar formation, 
respectively. The elapsed time in the calculation is shown in each panel. 



panels in Figure [3 respectively. The top panel of Figure [71 
shows that the disk gradually increases in size and surface 
density with time. In addition, the figure shows that a disk- 
like structure of size ^10 AU already exists before and im- 
mediately after protostar formation (solid line). This indi- 
cates that the remnant of the first core changes directly into 



the circumstellar disk. Thus, the formation and evolution 
processes of the circumstellar disk in weakly turbulent cloud 
cores are qualitatively the same as those in slowly and r igidly 
rotating cloud cores seen in previous works (see, e.g., iBata 
l201ll :lMachida et al.ll2010l ). In the bottom panel of Figure 
a positive radial velocity arises in the range of 50-100 AU 
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at ~ 5.0 X 10^ years after the second collapse. This outward 
flow in the disk is due to the m — 2 mode of spiral arms. 

Figure [8] shows the Fourier amplitudes of the surface 
density at 100 AU for 2 x 10^ years after protostar formation 
for model 6. For t < 5 x 10^ years, the disk radius is smaller 
than 100 AU. Thus, these amplitudes during this epoch orig- 
inate not from the circumstellar disk but from anisotropic 
accretion from the infalling envelope. After the circumstellar 
disk grows sufficiently, unlike the model which have strong 
turbulence, the m — 2 mode dominates the other modes, es- 
pecially for t > 1 X 10^ years, as shown in the bottom-right 
panel of Figure (6] 

As seen in pr evious studies (|Laughlin et al. 1 Il998l : 
iLodato Ric3l2005l ). without gas accretion onto the disk, 
the gravitational instability of the disk tends to develop 
the m — 2 spiral arms. Because of the weak turbulence 
in this model, the gas falls almost isotropically onto the 
circumstellar disk. Thus, the prominence of the m — 2 
mode is attributed to the gravitational instability of the 
disk. On the other hand , the model with the strong tur- 
bulence, the m — 1 density perturbation mode is domi- 
nant in the disk. This m — 1 mode is considered to be 
caused by anisotropic gas accretion onto the circumstel- 
lar disk. Note that the initial density perturbation of the 
disk st rongly affects the non- line ar development of the spiral 
arms (|Laughlin k, Bodenheimerlll994l l. These major modes 
play an important role in angular momentum transfer and 
readjust the surface density toward a more stable configu- 
ration. 



3.2 Evolution of Disk Orientation 

Figure [9] shows the time evolution of the angle of orien- 
tation of the disk angular momentum with respect to the 
angular momentum of the entire initial cloud core for non- 
fragmentation models (models 1-6). Because the disk is 
mainly supported by rotation, the orientation of the disk 
angular momentum roughly corresponds to the disk orien- 
tation. The disk angular momentum (disk orientation) is 
parallel to the angular momentum of the initial cloud core 
when ^ = 0. In the figure, the origin of time is set to the 
epoch when the protostar forms. This figure shows that the 
disk orientation is generally misaligned with respect to the 
total angular momentum of its host core at its formation 
and dynamically changes during the main accretion phase 
regardless of the strength of the turbulence. In each model, 
the angle of orientation of the disk is 30°-60° toward the an- 
gular momentum of the initial cloud core at the protostar- 
formation epoch (t = 0). As the accretion proceeds, the 
angle of orientation gradually decreases, reaching < 15° 
within several 10^ years for all models. It is expected that 
the disk orientation will be aligned with the angular mo- 
mentum of the cloud core when the gas accretion ceases. 

The right panel of Figure [9] shows that the angle of 
orientation for model 2 (a = 0.6 and 7turb = 0.06) varies 
irregularly with time. In the model, a very compact disk 
(< 10 AU) forms at first. It accretes onto the protostar and 
disappears in a short duration of ^ 10^ years. Then, an- 
other larger disk forms by the subsequent mass accretion. 
Reflecting this disk disappearance, model 2 shows an irreg- 
ular evolution of disk orientation. This irregular evolution 
in the early disk evolution stage may be unrealistic because 
of artificial treatment of the accretion onto the sink parti- 
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cle. However, this treatment for sink particle does not qual- 
itatively change our claim that the disk orientation angle 
dynamically changes with time. 

3.3 Fragmentation and Evolution of Disks in 
Binary System 

In models 7-9, fragmentation of filament occurs and a binary 
or multiple stellar system appears, whereas only a single pro- 
tost ar forms in models 1-6. The fragmentation process in a 
turbulent cloud core differs considerably from that in a cloud 
core with initial rigid rotation. Figure [TOl shows the density 
distribution when the protostars form for models 7-9. The 
figure indicates that the binary or multiple system forms by 
fragmentation of the filament. The figure also shows that the 
disks in a binary or multiple systems are not aligned with 
each other. For these fragmentation models, fragmentation 
tends to occur near the edge of the filament. This fragmen- 
tation process is possible to be explained by the focal point 
scenario suggested by Burkert & Hartmann (2004). 

To investigate further evolution of a wide binary sys- 
tem and disk evolutions, we focus on model 7 (a = 0.2 
and 7turb = 0.3). Figure [TTI shows the time evolution of the 
binary separation and its orbital energy for model 7. The 
orbital energy is given as 



where r and v are the relative distance and relative velocity, 
respectively. Note that the mass of protostar changes during 
the simulation. The right panel shows that the orbital 
energy steadily decreases, indicating that the binary 
system loses its orbital energy by tidal interaction 
between the disks and stars. In the left panel, we can 
see the separation of the binary also decreases. There- 
fore, it is expected that the observed separation of a binary 
or multiple system does not reflect the initial separation. In 
the right panel, sharp increases of the orbital energy 
are caused by the close encounter. During the early 
phase of binary formation, each disk is more massive 
than its host star and the stars are also accelerated 
by their surrounding disks. In equation (|8]), since the 
kinetic and potential energies of the disks are not in- 
cluded, this orbital energy is not conserved. There- 
fore, the orbital energy shows a large fluctuation at 
every close encounter. Note that it is considerably 
difficult to estimate the orbital energy of whole sys- 
tem which consider the stars and disks, because the 
disks exchange their mass and a part of the disk gas 
is striped off when close encounter occurs. However, 
we can roughly understand the evolution of the sep- 
aration and the orbital energy of the binary system 
for a long duration with Figure llll Note also that, 
the binary stars sometime have a positive orbital en- 
ergy, which seems to indicate that the binary system 
is not gravitationally bounded. However, the sign of 
orbital energy has no meaning since we ignored the 
disk component in equation ([8|). We confirmed that 
the whole system (disks and stars) are gravitation- 
ally bounded during this epoch. 

Figure [T2l shows an edge-on snapshot around the binary 
system at the end of the simulation ( t = 2.5 x 10^ years 



after the second protostar formation). Because the disk ori- 
entation around each protostar is mainly determined by the 
local velocity field at which the protostar forms, the disks 
are not aligned with each other. The misalignment is also 
seen in model 8. Thus, it seems that the misalignment be- 
tween disks is general feature of binary /multiple systems 
which form via filament fragmentation and initially have 
large separation (> 100 AU). 

The left panel of Figure [13] shows the evolution of the 
orientation angle of the disks from the angular momentum 
direction of the host core and the mutual inclination between 
the primary and secondary disks. In this study, we call the 
circumstellar disk around the primary protostar the primary 
disk and that around the secondary protostar the secondary 
disk. The angle of the primary disk orientation is ^ 90° in 
the very early phase of disk formation and it decreases with 
time (right panel). The sudden increases at t = 6 x 10^, 
1.6 X 10^, and 2.4 x 10"^ years are due to close encounters, 
as seen in Figure llll At each close encounter, the angle 
variation for the secondary disk is larger than that for the 
primary disk, which is because of the secondary disk is less 
massive than the primary disk. Thus, the secondary disk 
is strongly affected by the tidal interaction. As seen in the 
right panel of Figure 1131 the primary and secondary disks 
are inclined from the orbital plane by approximately 50° and 
30° at the end of the simulations. 



4 DISCUSSION 

4.1 Dynamical Change of Disk Orientation 

As shown in Figure [9l the disk orientation is generally mis- 
aligned with the angular momentum of its host cloud core 
and dynamically changes during the main accretion phase 
even in a weakly turbulent cloud. This is because the angu- 
lar momentum of the entire cloud core is determined mainly 
by the velocity field, whose wavelength is comparable to the 
cloud scale, A ^ i^cioud- On the other hand, the angular 
momentum of the disk is determined by the local velocity 
field, whose wavelength is much smaller than the cloud size, 
A <^ i^cioud- There is no phase correlation between these 
scales. This nature is independent of the strength of turbu- 
lence. The disk orientation changes due to the mass accre- 
tion which brings the angular momentum. The variation in 
the disk orientation during the main accretion phase can ex- 
plain the precessing outflow (e.g., Cunningham et al. 2 0^9h 
because the outflow is believed to be driven by the circum- 
stellar disk. 



4.2 Misaligned Disks in Binary Systems 

When large-scale fragmentation occurs and the binary sys- 
tem initially has a wide separation, disks tend to be mis- 
aligned because there is no correlation between the velocity 
fields in different regions. Thus, when a wide binary system 
forms in a turbulent cloud core, disks are expected to be 
misaligned with each other and also misaligned with the bi- 
nary orbital plane at their formation epoch. Our simulation 
showed that t he misalignment can mainta in at least for sev- 
eral 10^ years. iRoccatagliata et aP (|201lh recently observed 
a proto-binary system embedded in a common envelope and 
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Figure 9. Time evolution of orientation of disk angular momentum. The orientation angle ^ = corresponds to the orientation of 
the angular momentum of the initial cloud core. Left panel shows models having ol = 0.4, in which the solid, dashed, and dotted lines 
correspond to models with parameters of 7turb =0.3, 0.1, and 0.06, respectively. Right panel shows models having a = 0.6, in which the 
solid, dashed, and dotted lines correspond to models with 7turb = O-l? 0.06, and 0.03, respectively. 
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Figure 10. Logarithm of face-on surface density when a binary or multiple system appears in models 7 (top left, 7turb = 0-3), 8 (top 
right, 7turb = 0-1) cind 9 (bottom, 7turb = 0.06). Models 7- 9 have the same initial thermal energy of o; = 0.2. 



showed that the disks are highly misaligned (edge-on and 
face-on disks). They suggested that such a binary system is 
formed by fragmentation of two different parts of the col- 
lapsing molecular cloud core. Our study can naturally ex- 
plain such misaligned disks in a binary system. Misaligned 
disks are also obse rved even in a multiple stellar system 
(|Ratzka et al.ll2009l ). Our result also indicates that such a 



system can form via fragmentation of the collapsing cloud 
core in a turbulent environment, as seen in Figure [TOl 

We can also expect the disk orientation in a binary 
syst em from ob s ervati ons of molecular outflows or optical 
jets. IChen et~al\ (|2008l ) observed a low- mass protostellar bi- 
nary system and discovered two high- velocity bipolar molec- 
ular outflows that were nearly perpendicular to each other, 
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Figure 11. Time evolution of the separation between primary and secondary protostars (left) and the orbital energy (right) for model 
7. 




Figure 12. Surface density distribution at end of calculation for model 7 {a = 0.2 and 7turb = 0.3). Calculation stops 2.5 x 10^ years 
after the secondary protostar forms. 



showing a quadruple morphology. They concluded that the 
disks in a wide binary system are not necessarily co-aligned 
after fragmentation. A quadrupolar morphology and mis- 
alignedoutflows are ofte n observed in youn^ binary sys- 
tems (iMizuno et a l. 1990: Gueth et al.ll200ll : lLee et al.ll2002l : 
ICarrasco-Go ^^z et aimO^ . Because outflows or jets are 
expected to be driven by the disk, misaligned disks are ex- 
pected to show misaligned outflows. Therefore, observations 
of wide binary systems and our results support the idea that 
a binary /multiple system forms via dynamical fragmenta- 
tion of the filament in a turbulent molecular cloud core. 

4.3 Further Disk Evolution and Implications for 
Planet Formation in Binary System 

Tidal interaction between disks and protostars in binary sys- 
tem is essential for further evolution o f disk orientation in 
binary system. IPapaloizou k TerauemI (|l995l ) showed that 
perturbation of the secondary star can align t he primary 
disk toward the binary orbital plane. iBate et"al] (2000) also 
showed that the disk can be aligned with the binary orbital 
plane on the precession time-scale owing to tidal interaction. 
Thus, whether the disk is aligned with the orbital plane de- 
pends on the parameters of the binary system (especially 



on the orbital period). On the other hand, observations of 
main-sequence binary systems indicate that the stellar rota- 
tional equatorial planes (which may reflect the orientations 
of the disks) in wide binary systems (> 40 AU) are often mis- 
aligned with the binary orbital plane (|Halelll994l ). Thus, in 
wide binary system, disk misalignment is often maintained 
even after the main accretion phase. 

When disk misalignment remains after the main accre- 
tion phase, it is important for planet formation and the or- 
bital evolution of planets in a binary system. When planet 
formation in a binary system is investigated, it is convention- 
ally as sumed that each disk is aligned with the binary orbita l 
plane (jMarzari Schollll2000l : iTsukamoto Makinoll2007l l. 
However, when the disk and the binary orbital plane are mu- 
tually misaligned, the planet formation processes drastically 
change. By the Kozai effect (iKozail 1 19621 1, inward migration 
of planet esi mals occurs. Thu s, a compact dense planet esimal 
disk forms (jXie et al.1 120111). On the other hand, t he Kozai 
migration mechanism ("Pabrvckv Tremainel l2007l ) seems to 
be preferable for explaining the existen ce of close-in p lanets 
if disk misalignment generally occurs. iNarita et al.1 (|2010h 
recently found that the orbital plane of the planets is gener- 
ally not aligned with the orbital plane of the binary system. 
This also seems to be explained by the idea that the disks 
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Figure 13. Left panel shows time evolution of angle of disk orientation with respect to angular momentum of initial cloud core for 
primary (solid line) and secondary (dashed line) disks, and mutual inclination between the primary and secondary disks (dotted line). 
Right panel shows the time evolution of the inclination of the primary (solid line) and secondary (dashed line) disks toward the binary 
orbital plane. 



were misaligned with the binary orbital plane during the 
formation epoch. 



4.4 Comparison with Disk Evolution Simulations 
in Turbulent Cloud Core 

iGoodwin et al. 1 (|2004al lbh investigated the evolution of tur- 
bulent cloud cores. As we mentioned in ^2.2.2[ they con- 
sidered an ensemble of simulations with different random 
number seeds, and showed that fragmentation frequently 
occurs even with a moderate accretion rate that is real- 
ized in cloud cores with initially weak turbulence {a = 0.45 
and 7turb = 0.05). Their results are qualitatively the same 
as ours: fragmentation occurs in turbulent cloud cores and 
turbulence promotes fragment ation. Howev er, the fragmen- 
tation condition shown in iGoodwin et al.l |2004al lbl) seems 
to be quantitatively different from ours, because fragmen- 
tation only occurs with high accretion rate [a < 0.2 and 
7turb > 0.06) in our simulations. This difference may orig- 
inate in the difference of initial conditions. In their study, 
the initial cloud cores with a Plummer-like density profile 
has a total mass of M = 5.4M0. On the other hand, in our 
study, the initial cloud cores with a uniform density profile 
has a total mass of M = IMq. Thus, although it is difficult 
to qua ntitatively compare our results with Goodwin et al.. 
(|2004al lbh. the difference in initial conditions may some- 
what changes the fragmentation condition. It is also pos- 
sible that the different numerical resolutions may also affect 
the fragmentation condition. They resolved the cloud core 
with 25000 particles, whereas we resolved it with 520000 
particles. Mass resolutions in their study and in our study 
are mres = 2.1 x 10~^ and 1.9 x 10~^ M©, respectively, 
which implies that spatial resolution in our simulation is 
a pproximately fi ve times higher than in theirs. As discussed 
in lNelsoiil (|2006l l, there is a possibility that the lower spatial 
reso lution may enhance disk fragmentation. 

IWalch et al.1 (2010') investigated the evolution of turbu- 
lent cloud cores whose mass is Mcore = 6.1 Mq. They fixed 
the mean Mach number of cloud cores to be = 1 and 
changed the minimum and maximum turbulent wavelength 
and random seeds. In their study, the sum of the thermal 
and kinetic energies exceeds the gravitational energy, (i.e.. 



<^ + 7turb > 1). Thus, initial cloud cores are gravitation- 
ally unbound as a whole. The cloud begins to collapse after 
the energy dissipates by radiative cooling. By this approach, 
the consistency of the density distribution and the velocity 
field can be realized before the collapse begins. But it is 
difficult to control the accretion rate onto the disk with sim- 
ple parameters. Note that, in our study, the accretion rate 
can be controlled by the parameters a and 7turb because we 
adopted cloud cores that were initially gravitationally unsta- 
ble. They pointed out that initial differences of the angular 
momentum rarely affects the evolution of the cloud and cir- 
cumstellar disk. Rather, a filamentary structure appearing 
in turbulent cloud cores affects the evolution of the protostar 
and disk. Their results seem to be in agreement with ours. 
However, whereas a clear spiral arm often develops after 
disk formation in our simulation, most of the disks formed 
in their simulations does not develop the spiral arms. This 
difference may come from difference in the accretion rate. 

The evolution of low-mass very cold cores with a cloud 
mass of Mcore = 1-28 Mq was investigated in IWalch et al.l 
(|2012l ). They focused on the relationship between the evo- 
lution of the core and the maximum wavelength of the tur- 
bulence. They pointed out that the dynamical fragmenta- 
tion of filaments is a major process for binary formation in 
turbulent cloud cores. They demonstrated that binary for- 
mation via disk fragmentation is rare in the early phase of 
protostar formation which seems to be in agreement with 
our results. Their initial conditions are very cold and turbu- 
lence is very weak. Thus, disks formed in their simulations 
are very compact compare d with ours (see. Fig. [1] or [5] and 
Fig. 2 of IWalch et al.]|2012l ). In this study, we focused on the 
relationship between the evolution of the circumstellar disk 
and the turbulent energy of the cloud core. Furthermore, we 
showed how the disk orientation varies in time which is not 
mentioned in their study. 

4.5 Validity of Initial Conditions 

4.5.1 Inconsistency between the velocity and density field 
and possible effects on the results 

As described in ^ we imposed the turbulent velocity field 
(Gaussian random field) on the uniform density field. Thus, 
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they are not self-consistent when the simulation begins. 
Thus, it is expected that the velocity field do work to es- 
tablish self-consistent density field during the early phase 
of cloud evolution and the turbulence would decay. This 
decay could change the effective turbulent energy for the 
simulations. Another possible effect is that the density field 
could be still inconsistent even after the disk formation 
because our initial conditions are gravitationally bounded 
(a + 7turb < 1) and gravitational collapse immediately be- 
gins. In this subsection, we discuss whether these effects 
change our results. 

In figure [131 to investigate how much the turbulence de- 
cay, we show the time evolution of the velocity dispersion in 
the isothermal phase (maximum density of cloud is less than 
the critical density, pmax < pc ) for model 4 as an example. 
In the figure, the velocity dispersion slightly decreases due 
to the loss of turbulence energy for t < 60000 years. Then, 
it turns to increase at t ^ 60000 years, and continues to 
increase for tgtrsimGOOOO years. This increase is attributed 
to the infall motion of the cloud. This figure indicates that 
the decrease of turbulent energy before the cloud collapse is 
not so large. Thus, it is expected that this energy loss rarely 
affects the subsequent cloud evolution significantly. We con- 
firmed that the turbulent energy loss in other models is also 
small. 

Next, we discuss how the inconsistent density field af- 
fects the evolution of the system. In this study, we showed 
(i) anisotropic accretion enhances the low-order mode of spi- 
ral arms in the disk, (ii) the disk orientation dynamically 
changes during the main accretion phase, and (iii) the ori- 
entation of the disks is mutually misaligned and also mis- 
aligned with the orbital plane in binary system if the binary 
stars form in different region (or initial separation is suffi- 
ciently large). They are all originated from the anisotropic 
gas accretion onto the disk. The inconsistent density field 
would be smoother than a more realistic (or self-consistent) 
density field. Thus, if the self-consistent density field is 
adopted for the initial conditions, it is expected that the 
anisotropic accretion may be enhanced Thus, all the above 
features should be more stressed but not suppressed. There- 
fore, our results are not changed by the inconsistent density 
field. 

4.5.2 Validity of High Mach Number or Highly 
Gravitationally Unstable Initial Conditions 

In this study, to comprehensively understand the relation 
between turbulence and disk evolution, we adopted over a 
wide range of parameters and investigated both the disk for- 
mation around single star and in binary system. To investi- 
gate the disk evolution process around the binary system, we 
adopted a bit artificial initial conditions: initial cloud core 
of model 7 has large Mach number (A4 > 1) which is much 
larger than the typical observed value in cloud cores with 
masses of order 1 Msolar, and cores of models 8 and 9 are 
in highly gravitationally unstable state (a + 7turb = 0.3 and 
0.26). If our results and conclusions are closely related to the 
artificial aspects of these initial conditions, it is problematic. 

However, we would like to emphasize that our finding is 
not related to it. The main results we obtained from the sim- 
ulations of model 7, 8 and 9 are that the disk orientation of 
binary system generally misaligned each other and from or- 



bital plane if the initial separation of protostars is sufficiently 
large. This is because there is no correlation between the lo- 
cal velocity fields where the protostars forms. It is fundamen- 
tal property of turbulence that there is no correlation be- 
tween the local velocity fields around the different positions. 
Thus, such a misalignment may occur in the actual molecu- 
lar cloud core in which binary system forms. This should be 
confirmed by the observatio n and we t hink the observations 
we mentioned above (Hale 1993 : iRoccatagliata et al...2oTll ) 
seem to support our results. 
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